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Introduction 


The analysis presented here is part of a larger effort to develop an 
analytic technique capable of predicting the landing characteristics of pro- 
posed aircraft configurations in the early stages of design. In this first 
analysis, a linear pilot- aircraft closed loop model is evaluated using exper- 
imental data generated with the NT- 33 variable stability in-fligjit simulator. 
The pilot dynamics are modeled as inner and outer servo loop closures around 
aircraft pitch attitude, and altitude rate-of-change respectively. The 
landing flare maneuver is of particular Interest as recent experience with 
military and other highly augmented vehicles has shown this teisk to be rela- 
tively demanding, and potentially a critical design point. A unique feature 
of the pilot model used here is the incorporation of an internal model of 
the pilot’s desired flight path for the flare maneuver. 


Model Development 


Data from Ref. 1 suggests the landing flare maneuver can be modeled as 
a closed-loop tracking task, as pictured in the following sketch. 



Everything inside the dotted box in this model pertains to the pilot. 

The pilot model is partitioned into two parts. The left hand part represents 
the pilot’s internal model of the flare maneuver, and provides the basic 
error signal for tracking. This part of the model will be discussed more in 
following paragraphs. The ri^t hand part of the model Includes all the 
pilot dynamics which, as will be discussed below, consist of two gains and a 
time delay, plus possible first order lead or lag compensation. 

As can be seen from the sketch, this model only considers pitch axis 
dynamics and assumes linearized small perturbation equations. The following 
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variables are pertinent to the discussion. 


e 

- rad. 

Aircraft Pitch Attitude 

• 

H 

- m/sec. 

Altitude rate of change, or vertical veloc: 
of aircraft 

» 

Hid 

- m/sec. 

The pilot's desired vertical velocity at 
touch-down 

He 

- m/sec. 

Error In desired vertical velocity as per- 
ceived by pilot 

1/Ty 

- 1/sec. 

Flare model Inverse time constant 

%h 

- rad sec/m 

Pilot's Internal gain to convert h to an 
attitude command 

6c 

- rad 

Pilot's Internal estimate of attitude re- 
quired to correct 

0e 

- rad 

Error between desired and actual aircraft 
attitude; as perceived by pilot 


- m/rad. 

Pilot's Internal gain to convert 6 to an 
elevator (pitch) control motion ® 

T 

- sec. 

Pilot's delay time; approximates Internal 
piocesslng plus neuromuscular del«^s 


- sec. 

Pilot's optional lead compensation time 

Tl 

- sec. 

Pilot's optional lag condensation time 


- N 

Pilot's pitch control force 


As suggested In Ref. 1, the pilot's Internal flare model, mentioned 
above, can be represented as a, first order exponential. In other words, as 
a linear relationship between H and H, depicted In the following phase plane 



define the pilot's desired sink rate along the flare trajectory as H^, %diere 
He - 1/Tp H + 

and 1/Tp, the Inverse tine constant of the corresponding exponential decay 
Is: 

AH Htd - Hp 
1/Tp « — 

AH 0 - Hp , 
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Data from Ref. 2 will be used to test the model described above. For the 
experiment described In Ref. 2, characteristics representative of modem 
fighter aircraft with augmented dynamics were modeled In the varld>le stabil- 
ity NT-33. A series of ILS approaches was flown using a nominal glideslope 
of .044 rad (2.5 degrees); each approach was concluded with a visual flare 
and landing task. Tlie flare maneuver was initiated at a nominal altitude 
of 15.24m (50 feet) above ground. A sink rate of .762 m/sec. (2.5 ft/sec.) 
at touchdown is considered acceptable. The configuration designated as 2-7 
In Ref. 2 will be used as a specific example for further discussion here. 

For the visual flare, it is assumed the pilot has selected a desired touch- 
down point on the runway so the important speed is equivalent ground speed, 
which for this example was 61.73 m/sec. (120 knots or 202 ft/sec.) at flare 
initiation. Using the definitions above, these values give the ideal flare 
model as, 

H - -.11(H + 6.92) m/sec. 
with the exponential solution, 

H - 22. 16 (e*^^*^-. 31) n 

The airplane, actuator, augmentation and feel system dynamics for configura- 
tion 2-7 of Ref. 2 were examined and it was determined that a reduced order 
model using only the short period and command augmentation dynamics plus 
feel system gain would be adequate in the freqiiency range of .1 ^ j<*> Jl 20. 
rad/sec. 

The inal closed loop model was generated by first closing an inner loop 
for atti ;ude control, and then an outer loop for altitude rate. Ref. 3 sug- 
gests bandwidth and phase angle criteria for attitude control during landing 
of 

wii = 1.2 rad/sec, » it/2 

‘ U)g 

with a pure time delay of 0.3 sec. For this analysis a pilot time delay of 
0.33 sec was used since this value gave a convenient first order Fade' ap- 
proximate expression for linear frequency domain analysis. The pilot's time 
delay and any compensation Introduced were lumped in the attitude loop 
closure. Parameters were selected for the inner loop to approach the cri- 
teria mentioned above as closely as possible, while maintaining reasonable 
closed loop damping, for example, Ccl * 0*5. For the configuration selected, 
this rationale resulted in a fairly low-gain system, but it was felt this 
would be consistent with the pilot's desire to not overcontrol so close to 
the ground. Figure 1 summarizes the characteristics of the attitude closed 
loop with no pilot compensation considered. Continuing with the ter loop 
closure, and accounting for Initial conditions as suggested in Rei. ' . the 
final closed loop flare model was. 
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where the lower case letters are used to indicate perturbation values. The 
final closed loop transfer function was, 

h -38.61 (S - 6) 

' ■ m ■ - — ■ — — ■ I. p — ■ I. 1 I 1 I I ■ — — ^ ■ I ^ II ■ ■ II ■ III I — ■ • 

Hp (S + 1.164) (S + 6.005) [S + .771 + j 1.43lHS + 8.356 + j 8.611] 

The forcing function for this model, using the conditions defined above, was 
a step scaled to represent the altitude rate at flare initiation. 


Model Evaluation 


Figure 2 summarizes the response characteristics of the final closed 
loop system, and Figure 3 compares the landing trajectories of the piloted 
closed loop model with an Ideal e}q>onentlal curve. The trajectories In 
Figure 3 representing the piloted system with lead and lag were generated 
from models developed by the same procedure as described above. The pilot 
lead and lag were Introduced as cascade compensation within the pitch atti- 
tude loop. Figure 3 shows the pure gain and lead cospensated models both 
result In early flares and overshoot the ideal exponential touchdown time, 
leading to long gliding landings which pilots mlg^t describe as "floaters". 
The lag conpensated model, on the other hand, generally follows the expo- 
nential, but with some low frequency oscillation. An exanq>le flare trajec- 
tory from the Inflight experiment Is also shown on Figure 3. As the figure 
shows, the real fll^t trajectory does not follow the exponential or the 
analytic model path very closely except in the last 4 seconds before kOuch- 
down, where the flight trajectory is somewhat similar to the exponential. 
Figure 4 presents a second coDf>arlson of the flight path characteristics, 
using the H vs H phase plane. This presentation clearly shows the oscilla- 
tory nature of the lag compensated closed loop model. The figure also 
clearly shows the difference between the actual flight trajectory and the 
ideal exponential. 


Planned Work 


Further analysis of the Inflight data Is planned to Investigate differ- 
ences between the experimental and predicted flare trajectories. It Is felt 
the closed loop model Is a valid approach for the landing flare analysis be- 
cause of pilot comments and the nature of observed pilot control inputs 
during the flare. Figure 5 shews some examples of the pilot's longitudinal 
control activity during flare. The experlmei .al data will be examined to 
determine if significant time-varying or nonlinear characteristics are pre- 
sent which would account for the discrepancies with the linear, constant co- 
efficient model. 

Reference 4 Investigated STOL aircraft approach and landing using the 
optimal control pilot model. The flare maneuver was treated as a time optl- 
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nal problem in that study* resulting In time varying control. An example 
flare trajectory from the Reference 4 results is plotted in Figure 4. Al- 
though the trajectory is rotated because of the steeper approach angle used 
in the STOL study, the shape generally resenbles the NT- 33 actual landing 
trajectoiry. This observation lends credence to the plan to extend the pre- 
sent investigation to consider time-varying closed loop models. 
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ncURE 1. Clo««d Loop Attltuda ■•■pofwc with Tura Cain Pilot 
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FIGURE 2* Closed Loop Flare Response with Pure Gain Pilot 
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FIGURE 3. Closed Loop Flare Trajectories 
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FIGURE 5. Pilot Control Inputs in Flare 
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